Fluctuation relation and heterogeneous superdiffusion in glassy transport 
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Current fluctuations and related steady state fluctuation relation are investigated in simple coarse- 
grained lattice-gas analogs of a non-Newtonian fluid driven by a constant and uniform force field, 
in two regimes of small entropy production. Non-Gaussian current fluctuations and deviations from 
fluctuation relation are observed and related to the existence of growing amorphous correlations and 
heterogeneous anomalous diffusion regimes. 



Fluctuation theorem (FT) and nonequilibrium work 
relations are results of remarkable generality representing 
an important step towards the formulation of statistical 
mechanics far from equilibrium [U, 0] ■ FT states that 
the ratio of probabilities of observing an entropy produc- 
tion Wt over a long time interval r, to that of observing 
the opposite value, —Wt, is 



Tlr{-Wr) 



(1) 



Nevertheless, the class of systems obeying Eq. ([T]) is un- 
known - even for stochastic dynamics, where FT is most 
easily derived [1, Q ~ for it is not clear a priori when 
the asymptotic large-deviation/long-time regime under- 
stood in Eq. ([T]) is attainable, and whether it does gen- 
erally reflect most situations of physical interest. In 
fact, genuine deviations from Eq. ([1]) have been early 
observed and are now well established in various con- 
texts % iSi, [m, [ii d • Their general characteriza- 
tion - if they are accidental in nature or rather bring rel- 
evant informations about the stationary measure - how- 
ever, remains a widely open problem. 

In this paper, I show that deviations respecting the 
time reversal invariance of Eq. Wt —Wt, gener- 
ally occur in a large class of stochastic dynamics, and 
are a signature of glassy correlations. Their origin is 
traced back to the presence of heterogeneous anomalous 
diffusion regimes that extend possibly beyond any range 
of physically accessible values of t. Two distinct small 
entropy production limits will be considered: (i) vanish- 
ing forcing, i.e., near equilibrium, and (ii) far from equi- 
librium when currents become small at increasing drive, 
i.e. a negative resistance regime. Deviations in the for- 
mer limit are induced by a transient subdiffusion, and 
tend to decrease at long times. On the contrary, devia- 
tions in the latter limit increase with t and are due to 
a long-lived superdiffusion regime. Against our naive ex- 
pectation based on the behavior of relaxational glassy 
systems (where dynamics is typically sub-diffusive in the 
aging regime) , we find that mean-square current fluctua- 
tions grow super-diffusively at high-density, even though 
the average current becomes vanishingly small at increas- 
ing field. This particularly surprising behaviour is sug- 
gested to occur generally in driven systems dominated 
by steric hindrance and cage effect. Systems of exper- 
imental relevance include shear-thickening/jamming flu- 



ids, where the glassy correlations are intimately related to 
dynamic heterogeneity. The existence of non-monotonic 
responses is a peculiar feature of nonequilibrium steady 
state (NESS) and gives us the possibility of testing FT 
beyond the linear-response, and in situations which are 
of physical interest for driven soft matter. 

The models we study are emerging as a new paradigm 
for the interpretation of glassy phenomena [IJI- When 
driven into a NESS they show a nonmonotonic depen- 
dence of the relaxation time on the applied force, a fea- 
ture which is qualitatively similar to the behaviour of vis- 
cosity in sheared concentrated suspensions [l^. Two in- 
gredients characterize the dynamics: (i) the cage effect-a 
universal feature of glassy systems, which is implemented 
on a coarse-grained scale through a local kinetic con- 
straint, and (ii) a nonconservative force which consists 
of a uniform and constant drive, allowing for nonzero net 
current in the NESS. For simplicity, we consider two- 
dimensional (2D) square lattice systems in which the 
force is applied along a lattice axes, and the kinetic con- 
straint takes the following form: a randomly chosen par- 
ticle can move to a randomly chosen nearest neighbour 
site if i) the site is empty and ii) the particle has at most 
two nearby particles before and after the move. Further, 
if the move of a mobile particle is attempted in direction 
opposite to the field, the probability is e^'^^, otherwise 
it is 1. 

The interplay of the two above ingredients gives a 
rather rich dynamic behaviour [Tsj . At small forces, the 
cage around a given particle is just slightly distorted 
allowing easily for particle flow. At increasing forces, 
however, particles become generally more caged by their 
neighbors, as moves against the field direction are much 
less probable. Local rearrangements are thus more dif- 
ficult and transport more obstructed. As we shall see, 
this situation leads to non-Gaussian current fluctuations, 
growing heterogeneous spatio-temporal correlations and 
anomalous diffusion regimes. 

For locally reversible and irreducible Markov chains, 
such as the one introduced above, the proof of FT is 
straightforward ,5|. One considers the action 



W^r(W)-^l0g 



w{at,crt+i) 
w{at+i,(Jt)' 



(2) 



where w{a,<7') > are the transition probabilities for 
jumping from configuration a to a' . If the "border term" 



2 




FIG. 1: 2D constrained exclusion process driven by a small applied field E — 0.1, at particle density p — 0.8. The time r is 
measured in MCs and L is the linear system-size. Left: PDF of current fluctuations in normalized unit. Full line represents the 
Gaussian. Middle: Temporal evolution of relative fluctuations for several system sizes. Right: Logarithmic probability ratio 
log(nr(J)/nr( — J)) vs. entropy production Wt = EJt (/3 — 1). Full line, with slope 1, represents the prediction of Eq. ([l}. 



B = log [/is(cr)/Ms(o'o)]i (where fig is the stationary mea- 
sure), is subextensively small in r, then the generating 
function of the action verifies (e~'^^^) = ^e"^^"'*'-'^^ ), 
which is just another form of FT. In particular, when /is 
is flat over fixed density configurations, then B = and 
FT holds at any time r. This special case occurs in the 
asymmetric simple exclusion process (ASEP) to which 
our model reduces in the absence of constraints. In the 
presence of constraints, however the NESS measure is 
not trivial: while the density profile is flat, particles are 
statistically more clustered in the transverse direction at 
larger field. This nonequilibrium fluctuation-induced at- 
traction appears in the transverse pair-correlation func- 
tion and is a consequence of the more hindered longitudi- 
nal transport at increasing field. Nonetheless, while the 
transverse diffusion slows down (due to the stronger at- 
traction) , longitudinal diffusion is enhanced and becomes 
superdiffusive over a growing range of times, at increas- 
ing field. Neglecting border terms in this situation is not 
generally allowed even for large r, and that is the ul- 
timate reason of the observed deviations from Eq. ([1]). 
The importance of border terms is discussed thoroughly 
in Ref. [loj . We now explore their physical implications 
for driven systems with glassy dynamics. 

Fluctuation relation.- To test Eq. ([T]) we perform 
Monte Carlo simulations of the above driven stochastic 
dynamics. A system of linear size L is initialized with a 
uniform distribution of particles at fixed particle density 
p, and is let to reach the NESS. Observables of interest 
are evaluated over time intervals of duration r, along a 
trajectory of motion lasting 10^ — 10^ Monte Carlo sweeps 
(MCs), depending on the system size. In particular, we 
consider the particle current J^, that is the signed num- 
ber of jumps over r along the applied field, and compute 
its probability density function (PDF), UriJr)- Since 

w{a,a') — 0(constraint) min{l, e"^^ ''''}, local detailed 
balance holds irrespective of local kinetic constraints. For 
mobile particles one has w{a,a')/w{a' ,a) — 1,6=*=^^, de- 
pending on the relative direction of unit displacement 
and applied field {E ■ dr = 0, ±E). Whereas for immo- 



bile particle, w{a,a')/w{a' ,a) = 1, no matter the value 
of E ■ dr. The action functional ^ can then be easily 
identified as Wr = (3EJr , consistently with the standard 
definition of entropy production obtained from the time- 
dependent Gibbs entropy formula [H]. Thus, computing 
Ht-^Jt) allows for a direct check of FT. Notice, that the 
absence of both potential and kinetic energy terms makes 
the PDFs of current, work and heat exactly identical. 

We have first checked that in the standard ASEP 
current fluctuations are Gaussian distributed and that 
Eq. ^ is obeyed for time as small as r = 1 MCs, as 
expected because of flat measure. Similar behavior is 
found for constrained driven dynamics at small density. 
At moderately higher densities and small fields, the par- 
ticle motion becomes weakly correlated and small non- 
Gaussian tails appear in the PDF of current fluctuations, 
see Fig.[T](left panel). Since the average current J = (Jr) 
does not depend on r, and the PDFs for various values 
of T fall on the top of each other when plotted in normal- 
ized units, one might naively expect that the asymptotic 
large-deviation regime has been attained. Actually, cur- 
rent fluctuations, a^, still retain a dependence on r, and 
increasing the system size makes it harder to reach the 
asymptotic value. Fig. [T] (middle). The slow decay of cur- 
rent fluctuations leads to deviations from Eq. ^ which 
become smaller and smaller at increasing r, see Fig. [1] 
(right). Though such deviations are expected to disap- 
pear at longer r, in fact recovering FT maybe difficult at 
larger density, because the Ohmic regime shrinks fl^ . 

On approaching the negative resistance transport 
regime, which is a simple rheological analog of shear- 
thickening behavior |15j . something more interesting 
takes place. First, current fluctuations becomes strongly 
non-Gaussian and asymmetrically distributed, see Fig. [2] 
(left). Similar asymmetric PDFs have been observed in 
a wide range of systems [ll]. Second, although the av- 
erage current quickly attains its asymptotic value, the 
PDF keeps evolving with t, as shown by the behaviour 
of high-order moments in Fig. [2] (middle). Skewness and 
kurtosis generally display a nonmonotonous dependence 
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FIG. 2: Same as in Fig. 1, but in the negative resistance regime, with a large applied field E = 2.8 and L — 45. Left: The full 
hne is a generalised Gumbel function, with fitting parameter a = 11, see [13. Middle: Temporal evolution of PDF high-order 
moments. Notice that relative fiuctuations are rescaled by a factor 65. For systems of linear size L — 45, the time-scale for 
observing current-reversal events is at most of the order of 10^ MCs. Right: Deviations from Eq. ^ here increase with t. 



on T, which is suggestive of that found in the local corre- 
lations during the aging dynamics of glassy systems , 
and can be qualitatively understood as follows. On short 
time scales most particles are blocked and only a small 
fraction of them is able to move in the field direction, 
thus the PDF is left-skewed. Local rearrangements in- 
volving backward steps, however, are needed to sustain 
the system flow. Indeed, on intermediate time scales, 
current reversal events become more probable and, cor- 
respondingly, the PDF is less asymmetric. On longer 
time scales, most particles have moved and the probabil- 
ity of observing negative current is pretty small (although 
current reversal events keep happening on shorter time 
scales). The PDF thus changes to a right-skewed shape 
before eventually reaching the Gaussian form typical of 
diffusive, uncorrelated motion. Figure [2] (middle) shows 
that the latter diffusive regime occurs at times that are 
much beyond the time-scale required for observing neg- 
ative currents (at least two decades larger in the most 
favourable case, corresponding to parameters of Fig. [2]). 
Therefore, the asymptotic regime of large-deviation un- 
derstood in Eq. ([1]) is hardly attained, even for such small 
systems. What type of deviations should we expect in 
the present case and what features of NESS they pos- 
sibly encode? In spite of PDF asymmetry, we find that 
deviations from Eq. ([1]) are still well described by straight 
lines with a slope smaller than one, and their temporal 
evolution is mainly determined by current fluctuations, 
a^. However, unlike the previously discussed linear trans- 
port regime, the slope appears to decrease at longer t, 
see Fig. [5] (right). The numerical extrapolation of the 
asymptotic slope is difficult in this regime, because the 
accessible values of r are limited by the larger system-size 
used. In fact, sampling the PDF requires some care at 
high-density and large-force: to avoid finite-size induced 
dynamical blocking effects (leading to bimodal PDFs) 
one needs to use moderately larger systems. Since, this 
obviously decreases the probability of observing current 
reversal, we consider density barely larger than the neg- 
ative resistance threshold. Even in this most favorable 
situation (in which the saturation current is finite p^). 



current reversal events become unobservable well before 
the asymptotic regime is attained. The trend of high- 
order moments thus makes us confident that Eq. ([1]) will 
never be recovered in any realistic simulation. 

Heterogeneous anomalous diffusion.- Further insight 
into the nature of deviations from Eq. ([T]) is obtained by 
looking at the behavior of the longitudinal mean square 
displacement relative to the center of mass, Ar^ (t). Since 

the latter quantity is proportional to ct^, the physical ori- 
gin of the above behaviour can be traced back to the ex- 
istence of heterogeneous anomalous diffusion regimes, see 
Fig. [31 In the linear transport regime, Ar^(t)/t first de- 
creases with t, and then reaches an asymptotic diffusive 
plateau, that increases with the applied field. This en- 
hanced diffusion behaviour is also displayed by a driven 
particle probing an equilibrium glassy environment p^ . 
On approaching the negative resistance regime, the early 
subdiffusion range tends to shrink and connects to the 
late normal (enhanced) diffusion through an intermedi- 
ate (logarithmic) supcrdiffusive transient. Similar non- 
monotonic changes of the effective anomalous diffusion 
exponent have been observed in simulations of biased 
Brownian motion in a crowded environment [20], and 
are consistent with some continuous-time random-walk 
models [2l|. In the negative resistance regime, subdif- 
fusion is strongly suppressed, and one only observes a 
long-lived superdiffusion behavior, see Fig. [3] Although, 
we expect normal diffusion to occur eventually for any 
finite system, interestingly, we find that the crossover 
time to normal diffusion increases with the system size. 
The appearance of superdiffusion is particularly intrigu- 
ing in this context. It signals the simultaneous onset 
of long-time correlations in the motion of particles and 
the growth of their spatial correlations, as shown by 
the increasing peak of dynamical susceptibility [Tsj . A 
strikingly similar connection between anomalous diffu- 
sion and dynamic heterogeneity has been recently ob- 
served in vibrated granular packings, albeit in a dis- 
tinct jamming/rigidity regime [l^, see also [2^. It is 
tempting to conjecture, that the most peculiar feature 
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FIG. 3: Time averaged longitudinal mean-square particle dis- 
placement, Arj^|/f, vs. time t. 

distinguishing the approach to the near-jamming/shear- 
thickening/negative-resistance regime from the glassy be- 
havior, is that in the latter diffusivity tends to vanish, 
while in the former it grows unboundedly. 

Finally, it should be remarked that the results reported 
above are not specific of 2D: they were observed in 3D, for 
different direction of the applied force/lattice geometry, 
and for different types of kinetic constraints as well. 

To conclude, our results generally show that for a 
large class of driven stochastic dynamics, the asymptotic 
regime in which the steady-state FT holds, exceeds any 
reasonable physical time-scale, even for systems of mod- 



est size. The presence of heterogeneous superdifFusion 
in the negative resistance transport regime possibly sug- 
gests that the large-deviation function might not even 
exist in the usual sense, if the large-size limit is taken 
before the long-time one In spite of the above 

limitations, we find that deviations from ([1]) have sev- 
eral interesting features. First, they encode important 
physical properties, such as heterogeneous dynamics and 
anomalous diffusion. The strong "finite-time" deviations 
from ID) are a signature of growing long-range dynamical 
correlations both in time and space, and are arguably a 
universal feature inherent to driven systems with glassy 
dynamics, such as shear-thickening/jamming fluids. Sec- 
ond, their linear form is the simplest one respecting the 
time-reversal invariance of ([T]), though the asymptotic 
limit understood in FT is not representative of the phys- 
ical situation. Finally, they are highly suggestive of the 
notion of correlation-scale dependent effective tempera- 
ture [1^ . While the occurrence of longitudinal superdif- 
fusion seems to prevent the applicability of such an ap- 
pealing concept, we find that transverse fluctuation dy- 
namics is well described by a generalised Einstein relation 
(in some analogy to systems of driven vortices with ran- 
dom pin ning 12611 ) . Whether a modified form of Einstein 
relation [23, [2^ does hold for longitudinal fluctuations 
remains to be seen. 
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